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CN . We present analyses of the spatial and spectral evolution of hard X-ray emis- 

sion observed by RHESSI during the impulsive phase of an Ml. 7 flare on 2003 
November 13. In general, as expected, the loop top (LT) source dominates at low 
! energies while the footpoint (FP) sources dominate the high energy emission. At 
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intermediate energies, both the LT and FPs may be seen, but during certain in- 



tervals emission from the legs of the loop dominates, in contrast to the commonly 
observed LT and FP emission. The hard X-ray emission tends to rise above the 
FPs and eventually merge into a single LT source. This evolution starts first 
at low energies and proceeds to higher energies. The spectrum of the resultant 



LT source becomes more and more dominated by a thermal component with an 
increasing emission measure as the flare proceeds. The soft and hard X-rays 
Q-V show a Neupert-type behavior. With a non-thermal bremsstrahlung model the 

brightness profile along the loop is used to determine the density profile and its 
evolution, which reveals a gradual increase of the gas density in the loop. These 
results are evidence for chromospheric evaporation and are consistent with the 
qualitative features of hydrodynamic simulations of this phenomenon. However, 
some observed source morphology and its evolution cannot be accounted for by 
previous simulations. Therefore simulations with more realistic physical condi- 
tions are required to explain the results and the particle acceleration and plasma 
heating processes. 

Subject headings: acceleration of particles — Sun: chromosphere — Sun: flares — 
Sun: X-rays 
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1. Introduction 

Chromospheric evaporation was first suggested by Neupert (1968) to explain the origin 
of the hot, dense, soft X-ray (SXR) emitting plasma confined in the coronal loops during 
solar flares. The basic scenario is as follows. Magnetic reconnection, believed to be the 
primary energy release mechanism, heats the plasma and accelerates particles high in the 
corona. The released energy is transported downward along the newly reconnected closed 
flaring loop by non-thermal particles and/or thermal conduction, heating the chromospheric 
material rapidly (at a rate faster than the radiative and conductive cooling rates) up to a 
temperature of ~10 7 K. The resulting overpressure drives a mass flow upward along the loop 
at a speed of a few hundreds of km s -1 , which fills the flaring loop with a hot plasma giving 
rise to the gradual evolution of SXR emission. This process should also result in a derivative 
of the SXR light curve in its rising portion that closely matches the hard X-ray (HXR) light 
curve, which is called Neupert effect and is observed in some (but not all) flares (Neupert 
1968; Hudson 1991; Dennis & Zarro 1993; Dennis et al. 2003; Veronig et al. 2005). 

Hydrodynamic (HD) simulations of chromospheric evaporation have been carried out 
with an assumed energy transport mechanism (e.g., electron "beam" or conductive heating) 
(Fisher et al. 1985a; Mariska et al. 1989; Gan et al. 1995; Yokoyama & Shibata 2001; Allred 
et al. 2005) leading to various predictions on the UV-SXR spectral lines produced by the 
evaporated plasma, as well as the density and temperature profiles along the flaring loop. 
Most of the observational tests of these predictions rely on the blue-shifted components of 
SXR emission lines produced by the up-flowing plasma, first reported by Doschek et al. 
(1980) and Feldman et al. (1980) using spectra obtained from the P78-1 spacecraft. Similar 
observations were subsequently obtained from X-ray spectrometers on the Solar Maximum 
Mission (SMM; Antonucci et al. 1982, 1984), the Hinotori spacecraft (Watanabe 1990), the 
Yohkoh spacecraft (Wulser et al. 1994), and the Solar and Heliospheric Observatory (SoHO; 
Brosius 2003; Brosius & Philips 2004). Wulser et al. (1994), on the other hand, observed 
co-spatial SXR blue-shifts (up-flows) and H Q red-shifts (down-flows) as expected from HD 
simulations (Fisher et al. 1985b). A summary of relevant observations from SMM can be 
found in Antonucci et al. (1999). 

All the aforementioned observations, however, were indirect evidence in the sense that 
the evaporation process was not imaged directly. Based on HD simulations, Peres & Reale 
(1993) derived the expected X-ray brightness profile across the evaporation front and sug- 
gested that Yohkoh/ SXT or X-ray imagers with equivalent or better spatial and temporal 
resolution should be able to detect the front. Indeed, Silva et al. (1997) found that the HXR 
and SXR sources of the 1994 June 30 flare moved toward the loop top (LT) during the impul- 
sive phase. Since the flare was located near the center of the solar disk, they identified such 
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motions as the horizontal counterpart of the line-of-sight motion revealed by the blue-shifted 
emission lines observed simultaneously by Yohkoh/BCS. 

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) with its supe- 
rior spatial, temporal, and spectral resolution (Lin et al. 2002) provides us with opportunities 
to study the chromospheric evaporation process in unprecedented detail. We report in this 
paper our analyses of the spatial and spectral evolution of a simple flare on 2003 Novem- 
ber 13 with excellent RHESSI coverage. Because the flare occurred near the solar limb it 
presented minimum projection effects and a well-defined loop geometry that allows direct 
imaging of the HXR brightness profile along the loop. The observations and data analyses 
are presented in §2, followed by a derivation of the evolution of the density profile along the 
flaring loop in §3. We summarize the major findings of this paper and draw conclusions in 
§4. 



2. Observations and Data Analyses 

The flare under study is a GOES M1.7-class flare that occurred on 2003 November 13 
in the Active Region 0501 after it appeared on the east limb. This event followed a period of 
extremely high solar activities in late October and early November when a series of X-class 
flares including the record setting X28 flare of 2003 November 4 took place (Xu et al. 2004; 
Liu et al. 2004; Metcalf et al. 2005; Veronig et al. 2006). RHESSI had an excellent coverage 
of this flare. Figure 1 shows the RHESSI and GOES-10 light curves. The GOES 8-1 A 
(1.6-12.4 keV) and 4.0-0.5 A (3.1-24.8 keV) fluxes rise gradually and peak at 05:00:51 and 
05:00:15 UT, respectively. The RHESSI high energy (> 25 keV) count rates, on the other 
hand, exhibit two pulses peaking at 04:58:46 and 05:00:34 UT with the first one stronger. 
The steps in the RHESSI light curves are due to the attenuator (shutter) movements (Lin 
et al. 2002). Before 04:57:57 UT and after 05:08:59 UT, there were no attenuators in and 
between the two times the thin attenuator was in except for a short period near 05:05 UT 
when the attenuator briefly moved out. 

Figure 2 shows the evolution of the flare at different energies, which may be divided into 
three phases: (1) before 04:57:57 UT, corresponding to the rising phase when the emission 
mainly comes from a flaring loop to the south. (2) Between 04:57:57 and 05:08:59 UT, the 
impulsive phase, during which another loop to the north dominates the emission. This loop 
appears to share its southern footpoint (FP) with the loop to the south, which is barely visible 
because of its faintness as compared with the northern loop and RHESSIs limited dynamic 
range of ~10. (3) After 05:08:59 UT, the decay phase, when the shutters are out and two off- 
limb sources (identified as the LTs of the two loops) dominate. The relatively higher altitudes 
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compared with earlier LT positions are consequences of the preceding magnetic reconnection, 
as seen in several other RHESSI flares (Liu et al. 2004; Sui et al. 2004). Clearly the southern 
loop, which extends to a relatively higher altitude, evolves slower and is less energetic than 
the northern one. We shall focus on the evolution of the northern loop during the first HXR 
pulse (04:58-05:00 UT) in this paper. 

It is necessary to check if pulse pileup 1 is important in this flare, before we can make 
more quantitative interpretation of the data. The reason is that although we have applied 
the first order pileup correction (Smith et al. 2002) in our spectral analysis, such a correction 
is challenging for images and not available at present. There are several ways to do the check, 
the detector livetime being the first and simplest indicator. We first accumulated spatially 
integrated spectra for every 1 s time bins during the interval of 04:58:01-04:59:49 UT 2 , using 
the front segments of all the nine detectors except detectors 2 and 7 which have degraded 
energy resolution (Smith et al. 2002). We then obtained the livetime (between data gaps) 
from the spectrum object data and averaged it over the seven detectors being used. The 
resulting livetime generally decreases with time, ranging from 96% to 89%, with a small 
modulation produced by the spacecraft spin. In this Ml. 7 flare, such a livetime is comparably 
high (cf. livetime of ~ 55% during the 2002 July 23 X4.8 flare and ~ 94% during the 2002 
February 20 C7.5 flare) and indicates minor pileup severity. 

Another approach involves inspecting the change of the spectrum due to pileup. We 
accumulated spectra over every 1 spacecraft spin period (~ 4 s, with the same set of detectors 
mentioned above) and used the pileup correction to obtain the relative fraction of the pileup 
counts among the total counts as a function of energy (Smith et al. 2002). We find that 
the pileup counts amount to less than ~ 10% of the total counts at all the energies until 
04:59:01 UT when the livetime drops to 91%. After that, the relative importance of the pileup 
counts continues increasing, but remains below ~ 20% of the total counts before 04:59:17 UT. 
Toward the end of the first HXR pulse (04:59:45-04:59:49 UT, livetime ~ 90%), the pileup 
counts to total counts ratio exceeds 10% in the entire 20-40 keV range and humps up to 43% 
near 28 keV. We respectively integrate the pileup counts and total counts over the 20-40 keV 
band, and plot their ratio versus time, as a general indicator of pileup severity (see Figure 
3). Clearly this ratio is < 15% during the first 2/3 of the interval shown and does not reach 
the moderate ~ 25% level until the very end. 



1 Two photons close in time are detected as one photon, having their energies added. Pileup of three or 
more photons is possible but at a much lower probability (Smith et al. 2002). 

2 This time interval is also used in studying the evolution of the source morphology in §2.1 (see text about 
Figure 6), which covers the bulk of the first HXR pulse. 
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We therefore conclude that pileup effects are generally not very significant for this flare, 
especially during the first minute of the impulsive phase, because the count rate is not too 
high and the thin shutter is in at times of interest which further attenuates the count rate. 
It should be noted that the two piled up photons (that result in a single photon seen in the 
image) most probably originate from the same location on the Sun and pileup of photons 
across difference sources is relatively unimportant (G. Hurford, private communication). 
Therefore the source geometry would not be significantly affected by pileup, except that 
there could be a "ghost" of a low-energy source appearing in a high-energy image for very 
large (e.g. X-class) flares. However, the spectra of individual sources derived from images 
are distorted, which is relatively more significant at the LT than at the FPs. This is because 
ample low-energy photons dominate in population over high-energy photons, having the 
highest probability to produce pileup, and generally most of the low-energy photons are 
emitted by the LT source. 



2.1. Source Structure and Evolution 

We now examine the images in greater detail. The top left panel of Figure 4 shows 
RHESSI CLEAN (Hurford et al. 2002) images of the northern loop at 9-12, 12-18, and 28-43 
keV for 04:58:22-04:58:26 UT. (Although the 4-second integration time is rather short, the 
image quality is reliable with a well-defined source structure.) At 9-12 keV the LT dominates 
and the emission extends towards the two FPs, which dominate the emission at 28-43 keV 
and above with the northern FP (N-FP) much brighter than the southern one (S-FP). One 
of the most interesting features of the source structure is that emission from the legs of the 
loop dominates at the intermediate energy (12-18 keV). Similar structures are also observed 
for several other time intervals during the first HXR pulse (see discussions below). We find 
emission from the legs is a transient phenomenon at intermediate energies, because when 
integrating over a long period and/or a broad energy band, the LT and/or FP sources become 
dominant. To our knowledge no images like this have been reported before. We attribute 
this in part to the relatively short integration time and RHESSFs high energy resolution. 

To be compared with observations at other wavelengths, the same images at 9-12 and 
28-43 keV (solid contours) are shown with the SoHO/EYT, MDI magnetogram, and MDI 
white light maps in the other three panels of Figure 4, where the dashed contours depict the 
southern loop at 6-9 keV for 04:57:40-04:57:52 UT. The EIT image at 04:59:01 UT (upper 
right panel) shows emission at 195 A co-spatial with the SXR emission from the northern 
loop. The brightest 195 A emission, an indicator of the highest differential emission measure 
(thus the highest density) at ~ 1.3 x 10 6 K, appears to be close to the N-FP which is also the 
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strongest FP in HXRs 3 . The lower left panel displays the X-ray emission along with the post- 
flare (05:57 UT) MDI magnetogram. This clearly shows that the northern loop straddles 
across a polarity reversal with the brighter N-FP associated with a stronger magnetic field 4 . 
The southern loop (dashed contours) is associated with a even weaker magnetic field. Here 
we show the MDI magnetogram recorded one hour after the flare's impulsive phase because 
during a flare there are many uncertainties in the magnetic field measurement. The lower 
right panel shows the MDI continuum map at 12:47 UT (about 8 hours after the flare), 
suggesting that the flare occurred above the lower sunspot region (dark area). Note that 
during this interval the sunspot has moved westward for about 4° in heliographic longitude. 
We do not plot the MDI white light map at the time of the flare because then the sunspot 
was nearly on the limb and barely visible. 

Next we consider the evolution of the northern loop. We notice that, as shown in 
the four columns for 04:58:00-04:59:20 UT (boxed by the dotted line) in Figure 2, the FPs 
initially appear at all energies but later on dominate only in the high energy bands, while 
the LT is first evident at low energies and becomes more and more prominent at relatively 
higher energies as indicated by the dashed diagonal line. The emission from the LT also 
extends towards the legs at intermediate energies and in a given energy band the emission 
concentrates more and more at the LT with time. These are expected to be common features 
of flares with a single loop because of chromospheric evaporation that can increase the plasma 
density in the loop making the LT dominate at progressively higher energies. However, 
because the 20 second integration time is relatively long, these images do not uncover the 
details of the evaporation process. To remedy this we have carried out three different but 
complementary analyses of the images with higher temporal or energy resolution. 

1. To study the source morphology change over short time intervals, we model the loop 
geometry and study the evolution of the HXR brightness profile along the loop. We first 
made CLEAN images in two energy bands of 6-9 5 and 50-100 keV over the time interval of 



3 EIT 195 A passband images have a relatively narrow temperature response range with a characteristic 
temperature of 1.3 x 10 6 K (see Dere et al. 2000, Fig. 12), and emission intensity would be lower for both 
higher and lower temperatures. 

4 Note that since this flare occurred near the solar limb, the line-of-sight magnetogram measures mainly 
the horizontal (parallel to the solar surface) component of the magnetic field. The vertical component is 
more relevant here because flaring loops are usually perpendicular to the surface. However, it would be 
reasonable to assume that the vertical component scales with the horizontal one, and the polarity reversal 
line in the latitudinal direction is essentially not subject to the line-of-sight projection effect, as it seems very 
likely here. 

5 Since the thin attenuator was in at that time, counts below 10 keV are likely dominated by photons 
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04:58:12-04:58:53 UT that covers the plateau portion of the first HXR pulse. From these two 
images we obtained the centroids (indicated by the white crosses in Figure 5 a) of the sources 
identified as the LT (6-9 keV) and two FPs (50-100 keV), respectively. Assuming a semi- 
circular loop that connects the three centroids, we located the center of the circle marked 
by the plus sign in Figure 5 a. The grey scale in Figure 5 a was obtained by superposition 6 of 
30 images (six 8-second intervals from 04:58:08 to 04:58:56 UT in five energy bands: 9-12, 
12-15, 15-20, 20-30, and 30-50 keV) reconstructed with the PIXON algorithm (Metcalf et 
al. 1996; Hurford et al. 2002). Figures 5b and 5c respectively show the intensity profiles 
perpendicular to and along the loop (averaged over the respective orthogonal directions). 
The inner and outer circles (at r = 8'.'0 and 15'/3) in Figure 5a show the positions of the 50% 
of the maximum intensity in Figure 5b. However, for inferring the intensity profile along 
the loop we use radially integrated flux down to the 5% level. This enables us to include as 
much source flux as possible (with little contamination from the southern loop). We define 
the mean of the radii at the 5% level as the radius of the central arc of the loop (the white 
dot-dashed line in Figure 5 a). 

With the above procedure, one can study the evolution of the brightness profile along 
the loop at different energies. Figure 6 shows the results obtained from PIXON images with 
an integration time of 1 spacecraft spin period (~4 seconds) from 04:58:01 to 04:59:49 UT 
for three energy bands (20-30, 15-20, and 12-15 keV). Using a simple algorithm we determine 
the local maxima whose slopes on both sides exceed some threshold value and mark them 
by filled circles. We compare each profile with its counterpart obtained from the CLEAN 
image (with the same imaging parameters), and use the RMS of their difference to estimate 
the uncertainty as indicated by the error bar near the right-hand end of the corresponding 
profile. For each panel, the RMS difference of all the profiles, as a measure of the overall 
uncertainty, is shown as the error bar in the upper-right corner. This uncertainty is about 



whose real energy is about 10 keV higher than the detected energy. This is due to strong absorption of lower 
energy (< 10 keV) photons by the attenuator and escape of the germanium K-shell fluorescence photons 
that are produced by photoelectrical absorption of higher energy (10-20 keV) photons in the germanium 
detector (see Smith et al. 2002, §5.2). However, for the flare under study, the 6-9 keV image most likely 
reveals the real LT morphology, because there are ample thermal photons at lower energies originating from 
the LT source and photons at slightly higher energies seem to come from the same location. 

6 Bccause we are interested in determining the average loop geometry during the first pulse when the low 
energy X-ray flux has changed dramatically, using this approach to map the loop will ensure a relatively 
uniform brightness profile along the whole loop by assigning equal weights to images at different energies. 
On the other hand, if one simply integrates over the entire time range of 04:58:08-04:58:56 UT and energy 
band of 9-50 keV, the source morphology will be dominated by the LT source that emits most of the photons 
at a later time and relatively lower energies, which may not properly depict the loop geometry during the 
HXR pulse. 
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10% for the three energy bands; as expected, it increases slightly at higher energies which 
have lower counts. 

Figure Qa displays the profile at 20-30 keV, which as expected (see Figure 2) shows 
emission from the two FPs with fairly constant positions until the very last stage when the 
LT emission becomes dominant 7 . At this stage, the S-FP becomes undetectable and the N- 
FP has moved very close to the LT. At the lower energy (15-20 keV, Figure 66) the maxima 
tend to drift toward the LT gradually and eventually merge into a single LT source. At the 
even lower energy (Figure 6 c) this trend becomes even more pronounced and the drift starts 
earlier except here the shift is not monotonic and there seems to be a lot of fluctuations. We 
also repeated the same analysis at a higher cadence (every 1 second, ~ 4 second integration 
interval) with both PIXON and CLEAN algorithms. The evolution of the resulting profiles 
(although oversampled and thus not independent for neighboring profiles) appears to be in 
line with that shown here at a 4-second cadence obtained with PIXON. The general trends 
of these results indicate that high energy HXR producing electrons lose their energy and 
emit bremsstrahlung photons higher and higher up in the loop as the flare progresses. This 
can come about simply by a gradual increase of the density in the loop, presumably due to 
evaporation of chromospheric plasma. From the general drift of the maxima we obtain a time 
scale (~10's of seconds) and a velocity of a few hundred km s _1 consistent with the sound 
speed or speed of slow magnetosonic waves. As stated above at low energies we see some 
deviations from the general trend, some of which do not appear to be random fluctuations. 
If so and if we take one of the evident shorter time scale trends shown by the dashed line in 
Figure 6c we obtain a large velocity 8 (~10 3 km s -1 ) which is comparable to the Alfven or 
fast magnetosonic wave speed. This may indicate that another outcome of energy deposition 
by non-thermal particles is the excitation of such modes which then propagate from the FPs 
to the LT and might be responsible for the circularly polarized zebra pattern observed in the 
radio band (Chernov et al. 2005). This, however, is highly speculative because the spatial 
resolution (~7") is not sufficiently high for us to trust the shorter time scale variation. The 
longer time scale general trend, however, is a fairly robust result. 

2. Instead of examining the source structure with high time resolution, we can inves- 
tigate it with higher energy resolution at longer integration intervals as tradeoff for good 
count statistics and image quality. To this end, we have made PIXON images during three 



7 As noted earlier, pulse pileup in the 20-40 keV range becomes relatively important at this very late 
stage, which means that a fraction of the 20-30 keV photons seen in the image are actually piled up photons 
at lower energies. 

8 cf. Among the highest observed up- flow velocities in chromospheric evaporation are about 10 3 km s _1 
(Antonucci et al. 1990) and 800 km s _1 (Doschck et al. 1994), obtained from blue-shifted Fe XXV spectra. 
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consecutive 24 s intervals starting from 04:58:00 in 20 energy bins within the 6-100 keV 
range. Figure 7 shows a sample of these images at 04:58:24-04:58:48 UT. Figures 8a-8c show 
the X-ray emission profile along the loop at different energies for the three intervals 9 . As 
in Figures 6a-Qc the high energy emission is dominated by the FPs but there is a decrease 
of the separation of the FPs with decreasing energies and with time. Again at later stages 
the LT dominates and the profile becomes a single hump. The general trend again suggests 
an increase of the gas density in the loop. At lower energies (< 15 keV), the profile is 
more complicated presumably due to many physical processes (in addition to chromospheric 
evaporation), such as thermal conduction and transport of high energy particles, thermal 
and non-thermal bremsstrahlung, wave excitation and propagation, wave-particle coupling, 
and even particle acceleration, which may be involved. We believe that a unified treatment 
of acceleration and HD processes with physical conditions close to the flare is required for 
interpreting these results to uncover the details. 

To quantify this aspect of the source structure evolution, we divided the loop into two 
halves as shown by the boxes in Figure 7 and calculated their emission centroids. The 
resulting centroids at the three times together with the central arc of the model loop are 
plotted in Figure 9. As can be seen, for each time interval the centroids are distributed along 
the loop with those at higher energies being further away from the LT, and the entire pattern 
shifts toward the LT with time. Figure 10 shows the centroid positions of the northern half 
of the loop (where the source motions are more evident) along and perpendicular to the loop 
during the three intervals. This again shows that higher energy emission is farther away 
from the LT and the centroids shifted towards the LT with time, but similarly there are 
some complicated patterns at low and intermediated energies. All these are consistent with 
the general picture proposed above for the chromospheric evaporation process. 

3. To further quantify the source motions, we obtained the brightness- weighted stan- 
dard deviation or the 2nd moment of the profiles. In general the moment measures the 
compactness of the overall emission, but does not yield the sizes of individual sources whose 
measurement is still challenging for RHESSI (Schmahl & Hurford 2002). Hence our attention 
should be paid to the general trend of the moment rather than its absolute values, which 
may be subject to large uncertainties and thus less meaningful. The moments of the profiles 
resulting from CLEAN images (in three energy bands over 8 s intervals) are plotted in Figure 
116. There is a general decrease of the moment with the decline starting earlier at lower 
energies. Such a decrease is expected if the two FPs move closer to each other. However, 
caution is required here because a decrease of this quantity could also come about by other 



9 Note that pileup effects, as discussed earlier, are insignificant during this period of time (see Figure 3). 
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causes, say by an increasing dominance of the brightest source. We therefore checked the 
original images and the corresponding profiles when interpreting our results. To estimate 
the uncertainty of the moment, for each energy band we repeated the calculation with differ- 
ent integration time (panel c). The resulting moments remain essentially unchanged and as 
excepted the fluctuations of the moment decrease with increasing integration time. We also 
plot in panel c the moment (solid line) obtained from PIXON images with an integration 
time interval of two spin-period (~ 8 s), which basically agrees with its CLEAN counterpart 
in the general trend. The gradual 10 decrease of the moment is consistent with the motion of 
the centroids of sources up along the legs of the loop, which can take place by a continuous 
increase of the gas density in the loop due to evaporation. 



2.2. Spectral Analysis 

Spectral analysis can be used to study the evaporation process as well. With an isother- 
mal plus a power law model, we fitted the spatially integrated RHESSI spectra down to 6 
keV (Smith et al. 2002) for every 8 second interval during the impulsive phase. The emission 
measure (EM) and temperature of the isothermal component (asterisk symbols) are plot- 
ted in panels d and e of Figure 11, respectively. The EM rises almost monotonically with 
time from 0.6 to 14.2 xlO 49 cm -3 . This translates into an increase of the plasma density 
(n = a/ EM/V) by a factor of ~5 assuming a constant volume V . The temperature remains 
almost constant with a trend of slight decrease with time. The EM and temperature derived 
from the GOES data (plus symbols) are also shown for comparison. In general the GOES 
results are smoother and the temperature increases monotonically but remains below that 
of the RHESSI, consistent with previous results (Holman et al. 2003). This is expected be- 
cause RHESSI is more sensitive to higher temperatures than GOES. However, surprisingly, 
the GOES emission measure is also lower than that of RHESSI as opposed to what is the 
case more generally (see Holman et al. 2003). It is not clear whether or not this is due to 
a problem related to the RHESSI calibration at low energies. Nevertheless, the continuous 
increase of the £M's at comparable rates does suggest a gradual increase of the plasma 
density. 

The best fit parameters of the power law component with a low energy cutoff are plotted 
in Figure 11/. The power-law index 7 (plus symbols) is anti-correlated with the high energy 
light curves (see Figure 11a) and shows a soft- hard- soft behavior. It starts with 4.43 at 
04:58:02 UT, drops to 3.82 at the impulsive peak (04:58:26 UT), and rises up to 7.12 at 



On the other hand the jumps (if real) of the moment may suggest a transient phenomenon. 



- 11 - 



04:59:46 UT. The high indexes (> 5) may be an indicator of high-temperature thermal 
rather than non-thermal emission. Thus in what follows we limit our analysis to time up to 
04:59:20 UT. The low-energy cutoff (asterisks) of the power law is about 15 keV and is near 
the intersection of the isothermal (exponential) and power-law components. 



2.3. The Neupert Effect 

The Neupert effect is commonly quoted as a manifestation of chromospheric evaporation 
(Dennis & Zarro 1993) and a simple energy argument (e.g. Li et al. 1993) is often used to 
account for the relationship between SXR and HXR fluxes {Fsxr and Fhxr) ■ I n the thick- 
target flare model, the non-thermal Fhxr represents the instantaneous energy deposition 
rate (S e ) by the electron beam precipitating to the chromosphere, but the thermal Fsxr is 
proportional to the cumulative energy deposited, i.e., the time- integral of S e . It naturally 
follows that the temporal derivative of the SXR flux, F S xr, should be related to F HX r- 

The simplest test of the Neupert effect is usually carried out by plotting Fsxr and 
Fhxr in some energy band. There are many reasons why a simple linear relationship would 
not be the case here. The first and most important is that S e is related to Fhxr through 
the bremsstrahlung yield function Y (F HX r = £<Y) which is not a constant and depends 
on the spectrum of the electrons or HXRs (see e.g. Petrosian 1973). Here the most crucial 
factor is the low energy cutoff (Ei) of the non-thermal electrons, but the spectral index also 
plays some role. The total yield of all the bremsstrahlung photons produced by a power law 
spectrum of electrons with energies above E 1 (in units of 511 keV) is 




and the yield of the photons whose energies are greater than E x is 

where a = 1/137, In A = 20 is the Coulomb logarithm, and 5 is the spectral index of the 
power-law electron flux. As shown in Figure 11/ both the low energy cutoff and spectral in- 
dex of the non-thermal emission vary during the pulse, indicating variations in the electron 
spectrum and thus breaking the linearity of the SXR-HXR relationship. Other factors which 
can also produce further deviations are energy deposition by protons (and other ions), by 
conduction and other possible ways of dissipation of energy than simply heating and evap- 
orating the chromospheric plasma by non-thermal electrons. A detailed treatment of the 
problem requires solutions of the combined transport and HD equations, which is beyond 
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the scope of this paper. Veronig et al. (2005) by inclusion of some of these effects in an 
approximate way found that the expected relationship was mostly not present in several 
RHESSI flares. Finally one must include the fact that the chromospheric response of SXR 
emission will be delayed by tens of seconds depending on the sound travel time (and its 
variation) and other factors. 

The flare under study has shown no indication of gamma-ray line emission which means 
that the contribution of protons most probably is small. In the currently most favorable 
model where the electrons are accelerated stochastically by turbulence (see e.g. Petrosian 
& Liu 2004) the turbulence can suppress heat conduction during the impulsive phase and 
possibly also during the decay phase (Jiang et al. 2006). Because there does not appear to 
be large changes in the shape of the loop during the impulsive phase other energy dissipation 
processes such as cooling by expansion may also be negligible. Assuming these to be the case 
we have performed the Neupert effect test in two ways, the first being the common practice 
of examining the relation between Fsxr and Fhxr- We then examine the relation between 
£ e and Fsxr by taking into account the variation of the bremsstrahlung yield. 

I. The temporal derivatives of the fluxes of the two GOES channels are shown in 
the bottom panel of Figure 1. As evident, during the rising portion of the GOES fluxes 
both channels' derivatives indeed match the first pulse of the RHESSI HXR light curves 
(> 25 keV), but not during the second weaker pulse (where the 1-8 A derivative shows some 
instrumental artifacts). This may be due to the fact that the Neupert effect of the second 
pulse is overwhelmed by the cooling of the hot plasma produced during the first stronger 
pulse. Nevertheless, the SXR light curves (of both GOES and RHESSI) exhibit slightly 
slower decay rate than expected from the first pulse alone. This most likely is the signature 
of the energy input by the second pulse, which slows down the decay of the first pulse alone. 

We note in passing that the SXR light curves start rising several minutes prior to the 
onset of the HXR impulsive phase. This is an indication of preheating of the plasma before 
production of a significant number of suprathermal electrons. The 6-12 keV curve rises faster 
than the GOES curves at lower photon energies, which is consistent with the picture that 
the primary energy release by reconnection occurs high in the corona where the relatively 
hotter plasma is heated before significant acceleration of electrons (as suggested in Petrosian 
& Liu 2004), and before transport of energy (by accelerated electrons or conduction) down 
the flare loop to lower atmosphere where cooler plasmas are heated subsequently and pro- 
duce the GOES flux. On the other hand, the increase of the SXR flux at the beginning is 
dominated by the southern loop, which shows little evidence of chromospheric evaporation. 
The phenomenon therefore may be a unique feature of this flare. 

To quantify the SXR- HXR relationship, we cross-correlated the RHESSI 30-50 keV pho- 
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ton energy flux (F 30 _5o, Figure 12 a) and the derivative of the GOES low energy channel flux 
(F SX r, Figure 12c) in the SXR rising phase (04:58:00-04:59:51 UT). The resulting Spear- 
man rank-order correlation coefficient (see Figure 12/), an indicator of an either linear or 
nonlinear correlation, shows a single hump with a maximum value of 0.91 (corresponding 
to a significance of ~10~ 13 ) at a time lag of 12 s. This suggests a delay of Fsxr relative 
to F 30 _5o, which is expected given the finite hydrodynamic response time (on the order of 
the sound travel time of ~20 s for a loop size of ~10 9 cm and T ~ 10 7 K) required for 
redistribution of the deposited energy. Such a delay is evident in the numerical simulations 
of Li et al. (1993) who in addition found that the density enhancement contributes more to 
the total SXR emissivity than the temperature increase for longer duration (> 30 s) HXR 
bursts during the decay phase. In Figure 12c?, we plot the two quantities with the GOES 
derivative shifted backward by 12 s to compensate the lag of their correlation. A linear 
regression (dotted line) gives F 30 _ 5Q = (1.95 ± 0.15)Fsxr — (3.68 ± 0.48) with an adjusted 
coefficient of determination (or so-called R squared) E^ d - = 0.81 close to 1 suggesting a good 
linear correlation. 

2. We also carried out the same analysis for the electron energy power £ e , assuming 
a thick-target model of power-law electrons with a low-energy cutoff of E\ = 25 keV. We 
first obtained the energy flux of all the photons with energies greater than Ei, F El , from the 
30-50 keV photon energy flux -F 30 _5o: 



roo 771—7+2 

= J e J(E)EdE = F^ —l—— 2 , (3 ) 



n 00 7+2 

where J(E) cx E' 1 is the photon flux distribution at the Sun (photons keV -1 s _1 ) which is 
obtained from spectrum fitting (see §2.2) and assumed to extend to infinity in energy space. 
We then calculated the power of the electrons by 

£ e = F El /Y El , (4) 

where the bremsstrahlung yield Y El is given by equation (2) 11 . The resulting £ e is plotted 
versus time and GOES derivative in panels b and e of Figure 12 respectively. The dotted line 
in panel e shows a linear fit (R 2 adj =0.49) to the data: £ e = (0.65±0.11)F 5Xi? + (1.88±0.34). 
The corresponding Spearman correlation coefficient has a peak value of 0.78 (significance 
~10~ 8 ) at a time lag of 3 s (Figure 12/). As evident, £ e yields no better correlation with 
F S xr than -F 30 _ 50 , which is similar to the conclusion reached by Veronig et al. (2005). During 



n We used more accurate results from numerical integration of equation (29) in Petrosian (1973) rather 
than the approximate equation (2) here. However, one can still use equation (2) with a simple correction 
factor of 0.0728 x (S — 4) + 1 in the range 4 < 5 < 9 to achieve an accuracy of < 1%. 
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the HXR decay phase (after 04:59:20 UT) the spectrum becomes softer (7 > 5) and £ e 
decreases much slower than -F30-50 since the bremsstrahlung yield (equation [2]) decreases 
with the spectral index. As noted above for these high spectral indexes, the emission might 
be thermal rather than non-thermal. The inferred electron power is thus highly uncertain 
for these times. 

As stated earlier, the total energy of the non-thermal electrons is very sensitive to the 
low-energy cutoff E\ which is generally not well determined (cf. Sui et al. 2005). We thus 
set Ei as a free parameter and repeated the above calculation for different values of Ei 
(ranging from 15 to 28 keV). We find that, as expected, the temporal £ e -FsxR relationship 
highly depends on E\. For a small E\ (< 20 keV), £ e keeps rising till ~04:59:50 UT (near the 
bottom of the -F 30 _5 light curve), which makes the £ e -F S xn correlation completely disappear. 
On the other hand, for a large E x (> 20 keV) the correlation is generally good during the 
impulsive pulse (through 04:59:10 UT) and the larger E 1 the better the correlation. This 
is because the conversion factor E^ 1+2 / (30~ 7+2 — 50~ 7+2 ) in equation (3) is an increasing 
(decreasing) function of the photon spectral index 7 if E\ is sufficiently small (large). For a 
small Ei, for example, the photon energy flux F El may have a somewhat large value in the 
valley of the -F30-50 light curve when 7 is high. In addition, during this time interval the 
bremsstrahlung yield Ye 1 becomes small since 5 is large (see equation [2]) and consequently 
this may result in a very large £ e by equation (4). 

As to the magnitude of the energy flux of non-thermal electrons, Fisher et al. (1985a) 
in their HD simulations found that the dynamics of the flare loop plasma is very sensitive 
to its value. For a low energy flux (< 10 10 ergs cm~ 2 s _1 ), the up-flow velocity of the 
evaporating plasma is ~10's of km s _1 ; for a high energy flux (> 3 x 10 10 ergs cm~ 2 s _1 ), 
a maximum up-flow velocity of ~100's of km s _1 can be produced. For the flare under 
study, we estimate the area of the cross-section of the loop to be Ai oop < 1.6 x 10 18 cm 2 , 
where the upper limit corresponds to the loop width determined by the 5% level in Figure 
5b. We read the maximum electron power of £ e ,max = 9.8 x 10 28 ergs s -1 from Figure 126, 
which is then divided by 2Ai oop (assuming a filling factor of unity) to yield the corresponding 
electron energy flux: f e ,max ^ 3.1 x 10 10 ergs cm -2 s _1 . The source velocity on the order 
of a few hundreds of km s -1 estimated in §2.1 is consistent with that predicted by Fisher 
et al. (1985a). For comparison, we note that Milligan et al. (2006) from RHESSI data also 
obtained an energy flux of > 4 x 10 10 ergs cur 2 s _1 for an M2.2 flare during which an up- 
flow velocity of ~230 km s^ 1 was inferred from simultaneous co-spatial SoHO/CDS Doppler 
observations. 

In summary, the GOES SXR flux derivative F S xr exhibits a Neupert-type linear correla- 
tion with the RHESSI HXR flux F 30 _5o during the first HXR pulse. However, unexpectedly, 



the correlation between the electron power £ e and Fsxr is not well established based on 
the simple analysis presented here, which suggests that a full HD treatment is needed to 
investigate the chromospheric evaporation phenomenon (see discussions in §4). 



For the 1994 June 20 disk flare, Silva et al. (1997) interpreted the moving SXR sources 
as thermal emission from the hot (~ 30 — 50 MK) plasma evaporated from the chromosphere 
based on the good agreement of the emission measure of the blue-shifted component and that 
of the SXR from the FPs. For the limb flare under study here, Doppler shift measurements 
are not available. Meanwhile, a purely thermal scenario would have difficulties in explaining 
the systematic shift of the centroids towards the FPs with increasing energies up to ~70 keV 
as shown in Figure 10. A non-thermal scenario appears more appropriate. That is, the ap- 
parent HXR FP structure and motions can result from a decrease in the stopping distance of 
the non-thermal electrons with decreasing energy and/or increasing ambient plasma density 
caused by the chromospheric evaporation (as noted earlier in §2.1). One can therefore derive 
the density distribution along the loop from the corresponding X-ray emission distributions 
(e.g., Figure 8), without any pre-assumed density model (c.f. Aschwanden et al. 2002). This 
approach is described as follows. 

For a power-law X-ray spectrum produced by an injected power-law electron spectrum, 
Leach (1984) obtained a simple empirical relation (also see Petrosian & Donaghy 1999, §2) 
for the X-ray intensity I(r,k) per unit photon energy k (in units of 511 keV) and column 
depth r (in units of 1/ [47rr q In A] = 5 x 10 22 cm -2 for r = 2.8 x 10~ 13 cm and In A = 20): 



where 7 and 5 (= 7 + 0.7) are the photon and electron spectral indexes, respectively, A is 
a constant normalization factor, and dr = nds, where s is the distance measured from the 
injection site. This equation quantifies the dependence of the emission profile (or source 
morphology) on the electron spectral index and column depth. In general, when S decreases 
(spectrum hardening), the intensity at a given photon energy rises (drops) at large (small) 
r's and thus the emission centroid shifts to larger r's. This is expected because for a harder 
spectrum, there are relatively more high energy electrons that can penetrate to larger column 
depths and produce relatively more bremsstrahlung photons there. The opposite will happen 
when the spectrum becomes softer. During the impulsive peak showing a soft-hard-soft 
behavior (see §2.2), one would expect that the emission centroids shift first away from and 
then back toward the LT (if the density in the loop stays constant). Knowing the spectral 



3. Loop Density Derivation 
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index, the emission profile can therefore yield critical information about the density variation 
in both space and time. 

To compare the above empirical relation with observations, we first integrate 7(r, k) 
over an energy range [k\ , k 2 ] , 



and then integrate J(r; ki, k 2 ) over r to obtain the cumulative emission, 
Fir; k lM - [ At; kl ,h M r = £ [l - (l + r* + i) 

where we have chosen 

' 7 



(6) 



1-5/2' 



k^dk, (7) 



A ={L k ~" dk ) =w^- (8) 

so that F(t = oo; ki, k 2 ) = 1. Comparison of F(t; fci, A^) with the observed emission profiles 
gives the column depth t(s) whose derivative with respect to s then gives the density profile 
along the loop. 

Specifically for this flare, we assume that the non-thermal electrons are injected at the 
LT indicated by the middle vertical dotted line in Figure 8 and denote the profile to the 
right-hand side of this line (i.e., along the northern half of the loop) as J bs(s; fci, k 2 ), where 
[ki, k 2 ] is the energy band of the profile. The observed cumulative emission is then given by, 

PS P$max 

F obs (s; ki, k 2 ) = / J ob8 (s;ki,k 2 )ds/ / J obs (s;k 1 ,k 2 )ds, (9) 
Jo Jo 

where s max (corresponding to r = 00) is the maximum distance considered and F b s (s; k±, k 2 ) 
has been properly normalized. Then r = r(s; ki, k 2 ) can be obtained by inverting 

F(r;k 1 ,k 2 )=F obs (s;k 1 ,k 2 ), (10) 

where the integration over k in equation (7) can be calculated numerically. 

It should noted that, however, not all the profiles in Figure 8 are suitable for this 
calculation, because low energy emission is dominated by a thermal component especially 
in the LT region and at later times. We thus restrict ourselves to the energy ranges of 
12-72, 13-72, and 17-72 keV, respectively for the three 24 s intervals. The lower bound is 
the energy above which the power- law component dominates over the thermal component, 
determined from fits to the spatially integrated spectrum for each interval as shown in Figure 
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13. Within these energy ranges, separate leg or FP sources rather than a single LT source 
can be identified in the corresponding image, which is morphologically consistent with the 
non-thermal nature of emission assumed here. To further minimize the contamination of the 
thermal emission in our analysis, we have excluded the LT portion of the emission profile in 
excess of the lowest local minimum (if it exists) between the LT and leg (or FP) sources. An 
example of this exclusion is illustrated with the hatched region in Figure 8 c for the 19-21 keV 
profile. This was done by simply replacing the profile values between the LT and the local 
minimum positions with the value at the minimum. 

We calculated r(s; ki, A^) for every emission profile within the energy ranges mentioned 
above for the three intervals in Figure 8, with the photon indexes, 7 = 4.46, 3.97, and 4.23 
respectively. From the geometric mean of the column depths obtained at different energies, 
f, we derived the density profile n(s) = df(s)/ds for each time interval. The results are 
shown in Figure 14, where we shall bear in mind that attention should be paid to the overall 
trend rather than the details of the density profile and its variation, because the profile here 
only spans about 3 times the resolution (~ 7") and thus is smoothed, making neighboring 
points not independent. As can be seen, between the 1st and 2nd intervals, the density 
increases dramatically in the lower part of the loop, while the density near the LT remains 
essentially unchanged. The density enhancement then shifts to the LT from the 2nd to the 
3rd interval. This indicates a mass flow from the chromosphere to the LT. The density in 
the whole loop is about doubled over the three intervals, which is roughly consistent with 
the density change inferred from the emission measure 12 (see Figure lid). These results are 
again compatible with the chromospheric evaporation picture discussed in §2.1. 



4. Conclusion and Discussion 

We have presented in this paper a study of RHESSI images and spectra of the 2003 
November 13 Ml. 7 flare. RHESSFs superior capabilities reveal great details of the HXR 
source morphology at different energies and its evolution during the impulsive phase. The 
main findings of this paper are as follows. (1) The energy dependent source morphology in 
general shows a gradual shift of emission from the LT to the FPs with increasing energies. 
Over some short integration intervals emission from the loop legs may dominate at inter- 
mediate energies. (2) The emission centroids move toward the LT along the loop during 
the rising and plateau portions of the impulsive phase. This motion starts at low energies 



12 From 04:58:12 through 04:59:00 UT, the RHESSI {GOES) emission measure rises by a factor of 5.3 (2.3) 
which translate to an increase of the density by a factor of 2.3 (1.5), assuming a constant volume. 
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and proceeds to high energies. We estimate the mean velocity of the motion to be hundreds 
of km s -1 , which agrees with the prediction of the hydrodynamic simulations by Fisher et 
al. (1985a). There are also shorter time scale variations that imply much higher velocities 
(~10 3 km s _1 ) but we are not certain if they are real because of instrumental limitations. 
(3) Fits to the spatially integrated RHESSI spectra with a thermal plus a power law model 
reveal a continuous increase of the emission measure while the temperature does not change 
significantly. The GOES data show a similar trend of the EM but a gradual increase of the 
temperature. (4) The time derivative of the GOES SXR flux is correlated with the RHESSI 
HXR flux with a peak correlation coefficient of 0.91 at a delay of 12 s in agreement with 
the general trend expected from the Neupert effect. However, the correlation between the 
electron power and the GOES derivative is no better than the SXR- HXR correlation. (5) 
From the observed brightness profiles we derive the spatial and temporal variation of the 
plasma density in the loop, assuming a non-thermal thick-target bremsstrahlung model. We 
find a continuous increase of the density, starting first at the FPs and legs and then reaching 
to the LT. All these results fit into a picture of continuous chromospheric evaporation caused 
by the deposition of energy of electrons accelerated during the impulsive phase. 

Several of the new features of this event (such as the leg emission at intermediate 
energies) may be common to many solar flares. Expanding the sample of flares of this 
kind will be very helpful to understand the underlying physical processes. The new findings 
reside near the limit of RHESSrs current temporal, spatial, and spectral resolution. As 
advanced imaging spectroscopy capabilities are being developed and spatial resolution is 
being improved in the RHESSI software (Hurford et al. 2002), it will be critical to obtain 
the spatially resolved photon spectrum along the loop. This will yield incisive clues to the 
nature of the moving X-ray sources and relevant energy transport mechanisms and will be 
useful to check the reality of the short time scale variations. 

There are several important questions that need to be further addressed in future obser- 
vational and theoretical investigations: (1) What is the nature of the moving X-ray sources? 
Could they be characterized as thermal emission from the evaporated hot plasma or non- 
thermal emission from the precipitating electrons, or a mixture of both? Could they be 
related to MHD waves or evaporation fronts? (2) What are the roles of different heating 
agents of the chromosphere, i.e., electron beams, thermal conduction, and/or direct heating 
by turbulence or plasma waves during the impulsive phase? 

We have pointed out to some of the many physical processes that come into play in 
answering such questions. Here we describe possible directions for future theoretical studies. 
We have shown that a more physical based test of the Neupert effect between the electron 
power and SXR flux derivative does not reveal a better correlation than the usual HXR vs 
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SXR derivative correlation. Although the observed source velocity agrees with those of HD 
simulations, there are some features that current simulations have not addressed. To answer 
these questions requires an updated numerical calculation where one combines the model 
of particle acceleration and transport with the HD simulation of the atmospheric response 
to energy deposition to form a unified picture of solar flares. For example, one can use 
the output electron spectrum from the stochastic particle acceleration model (Hamilton & 
Petrosian 1992; Miller et al. 1997; Park et al. 1997; Petrosian & Liu 2004) as the input to 
the transport and HD codes, rather than simply assume a power-law electron spectrum as 
previous HD simulations. Such a study can shed light on the relative importance of particle 
beams and thermal conduction in evaporating chromospheric plasma and the roles that 
MHD waves may play in heating the flaring plasma, in particular, addressing our tentative 
observation of the fast source motion which suggests possible presence of MHD waves in the 
flare loop. A better understanding of their propagation, damping, and excitation mechanisms 
is necessary for uncovering the energy release process during flares. This is particularly true 
in the context of the stochastic particle acceleration model. 
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Fig. 1. — Top: RHESSI and GOES-10 light curves. The RHESSI count rates are averaged 
over every 4 seconds, with scaling factors of 1, 1/4, 1/12, and 1/50 for the energy bands 6-12, 
12-25, 25-50, and 50-100 keV, respectively. The sharp steps in the RHESSI light curves are 
due to attenuator state changes, and the sudden drop of the 6-12 keV count rate near 05:24 
UT results from the spacecraft eclipse. The GOES fluxes in the bandpass of 8-1 A (1.6-12.4 
keV) and 4.0-0.5 A (3.1-24.8 keV) are in a cadence of 3 seconds. Bottom: Time derivative 
of the GOES fluxes. Note that the periodic spikes of the low energy channel after 05:00:24 
UT are calibration artifacts. 
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Fig. 2. — Mosaic of CLEAN images at different energies (rows) and times (columns). Contour 
levels are set at 40, 60, and 80% of the maximum brightness of each image. The front 
segments of detectors 3-6 and 8 were used for reconstructing these images and the others 
presented in this paper, yielding a spatial resolution of ~7". We selected the integration 
intervals to avoid the times when the attenuator state changed. The large dotted box encloses 
the images during the first pulse of the impulsive phase, and within this time interval the 
dashed diagonal line separates the frames showing double sources or an extended source from 
those with a compact single LT source. 




Fig. 3. — The ratio of pileup counts to total counts, both integrated over the 20-40 keV 
range in time bins of 1 spacecraft spin. 
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Fig. 4.— Top left: RHESSI images for 04:58:22-04:58:26 UT during the first HXR pulse. 
The background is the image at 9-12 keV. The contour levels are at 75 and 90% for 9-12 
keV, 70 and 90% for 12-18 keV, and 50, 60, & 80% for 28-43 keV. Top right an EIT 195 
A image at 04:59:01 UT, showing co-spatial EUV emission in the northern HXR loop. The 
solid contours are the same as in the top left panel at 9-12 and 28-43 keV except that the 
contour levels are 50 and 80% for the latter. A 6-9 keV RHESSI image (same as the second 
panel in the first row of Figure 2) for 04:57:40-04:57:52 UT is plotted as dashed contours 
(at 50, 70, 90% levels) which depict the southern loop. The same set of contours is plotted 
in the two bottom panels as well. Bottom left: an MDI magnetogram at 05:57 UT. The 
line-of-sight magnetic field in the map ranges from -351 G (black: away from the observer) 
to 455 G (white) with the FPs near the strong magnetic field regions. Bottom right: an MDI 
continuum map at 12:47 UT showing the sunspots. The heliographic grid spacing is 2°. 
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Fig. 5. — (a) Synthesized image obtained by superimposing 30 8-second images between 
04:58:08 and 04:58:56 UT for 5 energy bands: 9-12, 12-15, 15-20, 20-30, and 30-50 keV. 
The three crosses mark the LT and two FPs identified as the emission centroids of the corre- 
sponding sources in the 04:58:12-04:58:53 UT images at 6-9 keV and 50-100 keV, respectively. 
The solid lines represent the semi-circular model loop with the center of the circles marked 
by the plus symbol. The white dot-dashed line is the central arc (see below) of this loop 
and the diamond indicates the start point of the distance in panel c. (b) Radial brightness 
profile averaged along the loop, obtained from the image shown in panel a. The distance 
is measured from the center of the circles. The horizontal dashed line marks the 50% level 
of the maximum and the crossings of this line with the profile define the radii of the two 
solid semi-circles in panel a. The 5% level is represented by the horizontal dotted line. The 
vertical dot-dashed line denotes the radial position of the central arc of the loop, (c) Same 
as panel b, but for surface brightness along the loop's central arc, averaged perpendicular to 
the loop. The three vertical dotted lines mark the corresponding positions of the cross signs 
in panel a. 
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Fig. 6. — (a) left: Evolution of the 20-30 keV brightness profile along the loop in a cadence 
of 4 seconds starting at 04:58:03 UT. Each profile is normalized to its own maximum and 
has an integration time of 1 spacecraft spin period (~4 seconds) whose central time is used 
to label the vertical axis. The filled circles mark the local maxima and the three vertical 
lines are the same as those in Figure 5 c. The error bar on each curve indicates an estimated 
uncertainty of the profile and the stand-alone error bar in the upper-right corner represents 
the overall uncertainty (13%) of all the profiles, (b) middle and (c) right: Same as panel 
a but for 15-20 and 12-15 keV, with an overall uncertainty of 12% and 10%, respectively. 
With the dashed straight line in panel c, we estimate the speed of the emission maximum at 
~10 3 km s _1 . Note the slightly different scales among the three panels for the profiles and 
their error bars. 
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Fig. 7. — PIXON images at 04:58:24-04:58:48 UT in different energy bands. The overlaid 
boxes were used to divide the loop into halves to calculate the corresponding centroids. 
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Fig. 8. — (a) left: Brightness profiles (obtained in the same way as Figure 6) at different 
energies for the time interval of 04:58:00-24 UT. The vertical axis is the average photon 
energy (in logarithmic scale) of the energy band for the profile. Representative energy bands 
(in keV) are labeled above the corresponding profiles. The vertical dotted lines are the same 
as in Figures 5 and 6. (b) middle and (c) right: Same as panel a but for 04:58:24-48 and 
04:58:48-59:12 UT, respectively. The error bars are the uncertainties of the corresponding 
profiles. The overall uncertainties, as indicated by the stand-alone error bar in the upper- 
right corner of each panel (note different scales, similar to Figure 6), are 14%, 13%, and 14%, 
respectively. The hatched region in panel c represents the LT emission (19-21 keV) being 
removed for deriving the density distribution in Figure 14 (see text). 
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Fig. 9. — Centroids of the northern and southern halves of the loop at different energies for 
the three 24-second time intervals (same as those in Figures 8a-8c). Energy increases from 
dark to light-grey symbols. The dot-dashed line marks the central arc of the model loop 
(same as in Figure 5a). 
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Fig. 10. — Positions of the northern centroids projected along (panel a) and perpendicular 
to (panel b, note the different scales) the central arc (the line in Figure 9) of the loop. The 
distance in panel a is calculated from the average LT position, as shown in Figure 5 a. 
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Fig. 11. — (a) RHESSI light curves (demodulated to remove artificial periodicity caused 
by the spacecraft spin), (b) Evolution of the standard deviation of the brightness profiles 
along the loop in three different energy bands obtained from CLEAN images, (c) Same as b 
but in the 15-20 keV band with different integration time intervals indicated in the legend. 
The solid curve denotes the result from PIXON images with an ~ 8 second integration 
time interval, (d) and (e) Evolution of the emission measure (10 49 cm -3 ) and temperature 
(MK), respectively, of the thermal component of the spatially integrated RHESSI spectrum 
obtained from fits to a thermal plus a power law model and from thermal fits to the GOES 
spectrum. The GOES emission measure is scaled by a factor of 10. (f) The evolution of the 
power-law index and the low-energy cutoff of the RHESSI power law component. 





Fig. 12. — (a) The 30-50 keV photon energy flux (.F30-50) at the Sun inferred from the 
RHESSI observation at 1 AU assuming isotropic emission. The two vertical dotted lines 
outlines the time interval (04:58:00-04:59:51 UT) used for the cross-correlation analysis (see 
below), (b) Power (£ e ) of the power- law electrons with a low energy cutoff of 25 keV inferred 
from the photon energy flux assuming a thick-target model, (c) Same as panel a but for 
the derivative (F S xr) of the GOES low energy channel (1-8 A) flux, (d) F 30 ~bo versus F S xr 
(shifted back in time by 12 s to account for its delay as revealed by the cross-correlation 
analysis; see panel f) within the interval of 04:58:00-04:59:51 UT. The grey scale of the plus 
signs (connected by the solid lines) from dark to light indicates the time sequence. The 
dotted line is the best linear fit to the data, (e) Same as panel d but for £ e and Fsxr which 
is shifted back by 3 s in time, (f) The Spearman rank-order correlation coefficient R of the 
photon energy flux (electron power) and Fsxr plotted as a function of time lag of the latter 
relative to the former. The dotted lines mark the peak values of R = 0.91 and 0.78 at a lag 
of 12 and 3 s, respectively. 
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Fig. 13. — The vF v spatially integrated spectra for the three 24 s time intervals as indicated 
in the legend. From the top to the bottom, the 2nd and 3rd spectra are shifted downward by 
2 and 4 decades, respectively. The broken lines are the thermal and power-law components 
of the fits to the data, and the solid lines are the sum of the two components. The thermal 
and power-law components intersect at about 12, 13, and 17 keV, respectively for the three 
intervals, above which the power-law component dominates. 
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Fig. 14. — Averaged density profiles along the loop inferred from the HXR brightness profiles 
during the three time intervals. 



